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Abstract. In the framework of spatially averaged inhomogeneous cosmologies in 
classical General Relativity, effective Einstein equations govern the regional and the 
global dynamics of averaged scalar variables of cosmological models. A particular 
solution may be characterized by a cosmic equation of state. In this paper it is pointed 
out that a globally static averaged dust model is conceivable without employing a 
compensating cosmological constant. Much in the spirit of Einstein's original model 
we discuss consequences for the global, but also for the regional properties of this 
cosmology. We then consider the wider class of globally stationary cosmologies that are 
conceivable in the presented framework. All these models arc based on exact solutions 
of the averaged Einstein equations and provide examples of cosmologies in an out— 
of equilibrium state, which we characterize by an information-theoretical measure. 
It is shown that such cosmologies preserve high-magnitude kincmatical fluctuations 
and so tend to maintain their global properties. The same is true for a A— driven 
cosmos in such a state despite of exponential expansion. We outline relations to 
inflationary scenarios, and put the Dark Energy problem into perspective. Here, it 
is argued, on the grounds of the discussed cosmologies, that a classical explanation of 
Dark Energy through backreaction effects is theoretically conceivable, if the matter- 
dominated Universe emerged from a non-perturbative state in the vicinity of the 
stationary solution. We also discuss a number of caveats that furnish strong counter 
arguments in the framework of structure formation in a perturbed Friedmannian model. 



PACS numbers: 04.20.-q, 04.20.-Cv, 04.40.-b, 95.30.-k, 98.80.-Es, 98.80.-Jk 
1. Introduction 

The standard model of cosmology idealizes spatial sections in terms of constant 
curvature hypersurfaces, the matter and energy distributions being spatially constant. 
Friedmann's and Lemaitre's solutions of Einstein's equations are currently employed 
to describe the dynamics of the Universe as a whole. Furthermore, the assumption is 
made that the spatially averaged inhomogeneous cosmos is described by a member of 
this family of solutions, a conjecture that can only be proved for Newtonian cosmologies 
|22] . |32j . but certainly corresponds to a restricted choice in General Relativity. 
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Friedmann-Lemaitre cosmologies describe a time-dependent, locally isotropic 
and, hence, on a simply-connected domain homogeneous and isotropic expansion or 
contraction, respectively. Denoting the scale-factor by a(t), which can be defined 
through the volume of an arbitrary comoving domain within the space sections, 
a := (V/V^ti)) 1 / 3 , normalized by the volume at some initial time £j, we obtain from 
Einstein's equations (restricted throughout this paper to a dust continuum) the well- 
known "acceleration law": 

a AtcGqh _ A 

a 3 3 ' 1 ' 

with the homogeneous restmass density qh, and the cosmological constant A . Its first 
integral yields a global expansion law, Friedmann's differential equation: 

a 2 8nGg H , k A g H (U) M 

3 + o? 3 " ; 6H ~ a 3 ~ V(U)a 3 ' [ ) 

with the total conserved restmass M enclosed within the arbitrary domain, and an 
integration constant k that is related to the spatially constant Ricci curvature Rh of 
the space sections by Rh = 6k/a 2 ; k := i£ff(tj)/6| 

Throughout this paper we shall call a particular solution of Friedmann's differential 
equation (J2J) a Friedmannian model or a Hubble flow. The only possible static 
Friedmannian model, a(t) =: ag, = const., which has motivated Einstein [HI] 
to introduce the (purely parametric) cosmological term A, follows from (J2J) with the 
assumption that a = in some finite time interval, and therefore also a = 0, and from 
Eqs. (HJ) and ©: 

AttGoe = A ; := ^f. — = const. ; (3) 
_^ _ STrGgg _ A 

a| 3 3 ■ 1 j 



— = A ; a E = /A _^^ , (5) 



Combining these two equations we obtain: 

— = A ; a E = 7== — 

and since ^ > 0, A has to be positive and hence also the curvature parameter k. 
Globally, we may introduce the total restmass in the Einstein cosmos Me = QeVe 
with its total volume V E '■= a 3 E Vi, and the global scalar curvature Re '■= 6k/a 2 E . Since 
the curvature is spatially constant and positive, the volume can be calculated in terms 
of spherical space. In space units (adopting now the normalization of the curvature 
parameter to k — +1), the Riemannian volume of the Einstein cosmos is Ve = 2^ 2 c 3 a 3 E , 

| In Friedmannian cosmology it is common practice to normalize the curvature parameter to k := 0, ±1 
by a suitable rescaling of the expansion factor a(t) (then acquiring the dimension of time for units with 
c = 1) . Although we formulated the equations for a dimensionless scale factor (to be consistent with 
the general equations to be discussed below), all numerical estimates in this section and especially in 
Appendix A will be done with the implicit understanding that k is normalized and a has dimension of 
time. 
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i.e. larger than the volume of a Euclidean sphere with the same radius 4ir/3c 3 a E . 
However, note that a mass-preserving smoothing of the Einstein radius into a Euclidean 
geometry yields a corresponding Euclidean volume 47r/3c 3 a E ir 3 that is larger than the 
total volume of the spherical space of Einstein's cosmos. Since qe is constant, we may 
rewrite the above equations as follows: 

AnGM E = AV E ; ^y + 8tcGM e + AV E = . (6) 

AVe may be interpreted as the total Dark Energy in this model. 

Einstein's model requires a non-vanishing and positive cosmological constant to 
"balance" the total restmass content of the Universe exactly (including the radiation 
density and pressure, which we here consider as being negligible in the matter-dominated 
era). Note that, by including presssure, there also exists a particular static model 
with A = and equation of state pn = — 1/3qh jUJ, p. 383). The beauty of a closed 
spherical space, as emphasized by Eddington (30], Ch.II), is accompanied by its definite 
predictions, e.g. for known restmass density we can determine its size as a strong 
boundary condition for any further studies. 

To illustrate this we calculate in Appendix A the Einstein radius by extrapolating 
the values of the cosmological parameters, as fitted to the Friedmannian model on 
the scale of our Hubble volume. Such estimates are rather naive as will become clear 
later, and we shall come back to this discussion within the more general setting of an 
inhomogeneous globally static cosmos. 

Soon after the time when the Einstein static model was suggested, the observed 
redshifts of galaxies together with their interpretation as Doppler velocities indicated 
that the space defined by the galaxies in our environment is expanding, which led to 
abandon Einstein's model. This was actually a hasty decision, based on the restricted 
view that the global model by Einstein was asked to describe any patch of the Universe, 
whatever small it was. These early discussions were based on the observational situation 
at the time which, following Hubble's assessment |33], indicates that the seemingly 
uniform distribution of galaxies may already represent a "fair sample" of the Universe. 
As a consequence, these discussions were meant in a global sense: the observed - 
according to contemporary standards very small - patch of the Universe was considered 
representative for the whole. Actually, Eddington jHI] already pointed out that the 
observed expansion might be a regional property of the Universe rather than a global 
one; he said "that it is possible that the recession of the spirals is not the expansion 
theoretically predicted; it might be some local peculiarity masking a much smaller genuine 
expansion; but the temptation to identify the observed and the predicted expansions is 
very strong". 

Before the Einstein cosmos was disregarded as a reasonable description of the 
Universe, there were many discussions following Einstein's in 1917: Dingle [2E] and 
Tolman [H2], [EH] have pointed out that the Einstein cosmos will soon develop into a 
highly irregular universe. This discussion is still referred to as the "instability of the 
Einstein cosmos" . Eddington [3T] argued that the Universe may evolve starting out from 
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the Einstein cosmos, however, due to its instability, will start expanding or contracting, 
respectively. As an alternative, Lemaitre [HU| , [SI] advocated the singular "Big Bang 
solutions" which expand until the matter density has dropped below the cosmological 
term which henceforth dominates, resulting in an accelerating phase thereafter (see 
also 7(J| §3C for a review of these discussions). The nowadays favoured concordance 
model (cf. Appendix A), featuring a positive cosmological constant, describes such an 
evolution. We shall see that the ideas advanced at that time apply to the picture 
developed below, but their interpretation will be very different. In contrast to the 
historical flaw of realizing the instability of the Einstein cosmos, we nowadays view 
such instabilities as the origin of large-scale structure. All homogeneous world models 
are unstable including the Einstein cosmos. As we shall discuss in detail below, the 
global instability argument in the above form (i.e., within the class of homogeneous- 
isotropic models) does not apply, when the cosmology acquires the status of describing 
the average dynamics on the largest scales. 

Here, a disclaimer is in order: to advance a globally static, but regionally fluctuating 
cosmos as a viable model that could explain current observational results is premature. 
Instead, we revisit the ideas which led to the Einstein cosmos (the introduction of the 
cosmological constant) in light of a new framework and on the grounds of an ongoing 
discussion of the possibility that Dark Energy may be explained through "backreaction 
effects" of structure formation. Globally expanding, stationary cosmologies are also 
conceivable in this framework. It was pointed out in [TS| that the question whether the 
magnitude of backreaction effects is sufficient to explain acceleration of the observable 
Universe, and the question of whether a globally static or stationary cosmos bears 
physical justifications beyond a mere mathematical possibility are synonyms. 

A thorough physical investigation of fluctuation-supported static and stationary 
cosmologies must be based on more general matter models (as a next step perfect fluids 
and scalar fields) that allow to study dynamical scenarios of inflation and their exit 
details including the effect of radiation pressure. The corresponding effective equations 
are given in ^H] • In this respect the present investigation based on the effective equations 
for a dust continuum [T3] provides a useful showcase for the presentation of the basic 
ideas. The present investigation of particular exact solutions offers more insight into 
the general formalism of averaged inhomogeneous cosmologies, and at the same time 
proposes new families of cosmologies that enjoy significantly more freedom than a rigid 
Friedmannian cosmology. Furthermore, by applying these ideas to inflationary scenarios, 
we can understand the relevance of matter and curvature fluctuations for the description 
of the Early Universe, in particular the importance of the role played by a non-vanishing 
averaged scalar curvature. 
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2. Effective equations for inhomogeneous universe models 

2.1. Averaged equations 

For the sake of transparency we shall restrict all considerations in this paper to an 
irrotational dust continuum and recall a set of effective equations provided in ^H] • The 
ideas presented can be carried over to studies of inhomogeneous cosmologies covering 
the Early Universe and radiation-dominated epochs with the help of the more general 
effective equations developed in [Tflj . 

Given a foliation of spacetime into flow-orthogonal hypersurfaces (which is possible 
for irrotational dust) with the 3-metric gij in the line-element ds 2 = —dt 2 + gij dX % dX\ 
spatial averaging of any scalar field \I/ is a covariant operation and is defined by the 
simple averager: 

:=i / Jd 3 X , (7) 

» 25 JV 

with J := ^/det(<7y); is the metric of the spatial hypersurfaces, and X 1 are 
coordinates that are constant along flow lines, which are here spacetime geodesies. 
Following jT^j we define an effective scale factor by the volume of a simply-connected 
domain T> in a t-hyper surf ace, normalized by the volume of the initial domain Pj, 

We recall the fact that, for a restmass preserving domain T>, volume averaging of a 
scalar function does not commute with its time-evolution: 

(dt9) v -dt{9) v = (9) V (9) V -(W) V , (9) 

where 9 denotes the rate of expansion. Setting ^ = g we obtain a regional continuity 
equation reflecting the conservation of the total restmass within T>: 

d t M v = & d t (g) v + (0) v (g) v = . (10) 

Setting = 9 we can derive an effective equation for the spatially averaged expansion 
of the model, 

(9)v = ^ = = : 3H V , (11) 

where we defined an effective Hubble functional on X> (an overdot denotes partial time- 
derivative): inserting Raychaudhuri's evolution equation, 9 = A — 47rG^ — |# 2 — 2cr 2 
(with the rate of shear a 2 = hcr^cr^), into (jSJ) and using the effective scale-factor a© we 
obtain: 

3^ + 4nG^-A = Q v . (12) 

The first integral of the above equation is directly given by averaging the Hamiltonian 
constraint: 

avV _ 8nGMv_ _ A = _Qv , . 

a v J 3 Viof, 6 3 6' 1 ' 
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where the total restmass Mp, the averaged spatial Ricci scalar (1Z) V and the kinematical 
backreaction term Qv are domain-dependent and, except the mass, time-dependent 
functions. The backreaction source term is given by 

Qv := 2 (IJ) V - ~ (1)1 = \{{9~ (6) v f) v - 2 (a 2 ) v ; (14) 

here, I = Q t i and II — |[ (0*j) 2 — 0^0^] denote the principal scalar invariants of the 
expansion tensor, defined as minus the extrinsic curvature tensor := 0jj. In the 
second equality above it was split into kinematical invariants through the decomposition 
0y = \gij9 + with the rate of expansion 9 = Q\, and the shear tensor ov,-. (Note 
that vorticity is absent in the present model; we adopt the summation convention.) 

The time-derivative of the averaged Hamiltonian constraint ("L 7 ")) agrees with the 
averaged Raychaudhuri equation ("l"2|) by virtue of the following integrability condition: 



d t Q v + Q—Qv + d t (R) v + 2^ [K) v = , (15) 
av dv 

which we may write in the more compact form: 

\d t (Q v a%) + ^((44) = 0- (16) 
Formally integrating this condition yields: 

% - A f d t' Qv ^aUt') = I ( (K) v + Q v ) , (17) 

i.e., besides the total material mass Mp we have a further integral of motion given by 
the domain-dependent integration constant kx> that we may also write as follows: 



/"* d 
Qk v = Y v + 2 dt' Q v ^4(0 " Qva 2 v . 



In the above equation we introduced the functional Yd, which is a special case of 
the Yamabe functional, here Yd := {I^d^v^ m t nree dimensions (see, e.g., the case 
n = 3; <p = ip = 1 in [4J, p. 150), which itself is an integral of motion for vanishing Qd§. 
Eq. (|16j). having no Newtonian analogue, shows that the averaged intrinsic curvature 
and the averaged extrinsic curvature (encoded in the backreaction term) are dynamically 
coupled. Stating this genuinly relativistic property, we also note the surprising fact that, 
inserting (|17|) into (|13|) results in an equation that is formally equivalent to its Newtonian 
counterpart [2*2"] : 

d% + k v 8ttG (q) v A 1 f \ , d 2 



/ dt' Qp-^4(t') . (19) 



a v 3 3 3a|> j ti 

The effective scale-factor obeys the same equation as in Newtonian theory similar to the 

situation known for the homogeneous-isotropic case. Note that these effective equations 

also cover anisotropic inhomogeneous cosmologies ("221, Appendix B). 

§ The vanishing of Q-p on every scale is necessary and sufficient for a-p to be a (global) solution of 
Friedmann's differential equation. Therefore, a non-vanishing Qx> plays a key-role and justifies the 
name "kinematical backreaction". Y-p is also preserved in the special case where Qv oc a^ 6 , which is 
briefly discussed at the end of this section. 
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2.2. The cosmic quartet 

For the purpose of comparing the model variables with observations it is comfortable to 
introduce dimensionless average characteristics as follows (contrary to "H"] we use the 
notation for the curvature functional): 



v _ 8nGM v v _ _A_ v _ (K) v v _ Q v 

m -~w v .alHi ' A ""3if2 ' n '~ ml ' s '~ ml 



where we have employed the effective Hubble-functional H v (fTTj) that reduces to 
Hubble's function in the homogeneous-isotropic case. With these definitions the 
Hamiltonian constraint is written in the iconized form of a cosmic quartet: 

o£ + n£ + ng + nS = 1 ■ (21) 

These functionals, being scale-dependent, are dynamically related in a complex way||, 
unlike the situation in a Friedmannian model that features a global cosmic triangle 
p"] where the cosmological parameters interact trivially, e.g. for an Einstein-de Sitter 
cosmology (vanishing curvature parameter and vanishing A) the curvature parameter 
remains zero throughout the entire evolution. 

As in Friedmannian cosmology we may also define other functionals, as for example 
an effective deceleration parameter. 



&:=-^4 = in£-n? + 2nS > ( 22 ) 



as well as corresponding effective parameters for the third derivative of a© like the 
recently introduced state finders P (see also and references therein). 

We also note the useful general evolution equation (a slightly different version is 
given in ^H], Appendix B): 

4-$(n24) + \d t (nlal) -afl^ng + ng) (5^ + ^ + 2^) =0.(23) 

A clarifying remark concerning the definition of the "curvature parameter" is in 
order. In the corresponding Newtonian problem the following curvature and kinematical 
backreaction functionals have been introduced [2*3*] : 

rt a 

(24) 

In view of Eq. (jl9j) we may use the same functionals in place of those introduced 
above also in General Relativity. However, the physical averaged curvature is not 
associated with Qf. The functionals (f2*4*|) are related to the previously introduced ones 

by: nv + nlx = n% + n». 

|| Here, the detailed investigation of these functionals in the framework of Newtonian cosmology in [23] 
is useful to understand this remark. 
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2. 3. The cosmic equation of state 

The above equations can formally be recast into standard zero-curvature Friedmann 
equations with new effective sources [16j^[: 

& = ^-xL Q »-iL {n) » ; * = -ie^ <fo + 4^<^ • (25) 

3— = A - 47rG(fiS + 3p£) ; 3fl*=A + 87rG<& ; & + 3flj, (<& = . (26) 
a© 

Eqs. correspond to the equations (JI2J), (|T3*j) and (JTOJ) , respectively. In these 

equations we have translated all what has been said before into a Friedmannian setting 
and a specific form of the fluctuating sources. Note that the kinematical backreaction 
term Qx> itself obeys a stiff equation of state mimicking a dilatonic field in the fluid 
analogy (for further implications see jTB] and Subsect. 13. 6|) . 

Given an equation of state of the form pf s = j3 (g^,ax>) that relates the effective 
sources with a possible explicit dependence on the effective scale factor, the effective 
Friedmann equations can be solved (one of the equations (j2T?j) is redundant). 
Therefore, any question posed that is related to the evolution of scalar characteristics 
of inhomogeneous universe models may be "reduced" to finding the cosmic state on a 
given spatial scale. (Note, however, that an equation of state must not exist in general.) 
Although formally similar to the situation in Friedmannian cosmology, here the equation 
of state is dynamical and depends on details of the evolution of inhomogeneities. In 
general it describes non-equilibrium states. 

An example may illustrate the cosmic equation of state: we look at a particular 
exact solution of the averaged Einstein equations that was given in ( |15j Appendix B). 
From the integrability condition (fTH|) we directly infer that the pair of solutions, 

Qv = 9fL , {n)v = <M , (27) 

a v a v 
provides a special solution for which the averaged scalar curvature and the kinematical 
backreaction term decouple and evolve independently. The regional cosmic equation of 
state corresponding to this solution can be derived from (|23|) : 

^ - - 1 ~ \^ ( 9R ) 

with 7x := Q^(ti) / Q^(ti) and 72 := f^(ii)/f2g(£i). Asymptotically, for an expanding 
domain (\a v \ large), wf s tends to — |, a property that is also shared by a non-flat 
Friedmannian domain for which we obtain t( ; Friedmann — — 1/(1 + 73/a©) with 73 : = 
Vt^ n {ti) /VL^iti). These Friedmannian models are subcases of the above inhomogeneous 
solution (flg(ti) ^ 0; fi%(U) = 0); w Friedmann follows from (jSHJ by first multiplying with 

Note that in this representation of the effective equations p e ff just denotes a formal "pressure" : in 
the perfect fluid case with an inhomogeneous pressure function the foliation has to be differently chosen 
(the in general inhomogeneous lapse function is set equal to 1 here), and there is a further averaged 
pressure gradient term |16j . 
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&q{U) 7^ (it is understood that the scale-factor and all other quantities then no longer 
depend on V). 

We end this section with a note on the definition of the effective sources in 
Eq. f)25j) : there is some ambiguity in defining them. Firstly, it may sometimes be 
useful to incorporate A into the effective sources by defining gf SA := gf s — A/8nG and 
PeSA := Pes + A/87rG. Secondly, we might add the "constant curvature term" 3/cx>/a|, to 
the left -hand-side of the second equation in (|26j) : if we wish to do so, then the effective 
sources can be represented solely through the kinematical backreaction term Qx> and its 
time-integral. For this we have to exploit the "Newtonian form", Eq. (|19p. and would 
have to define the effective sources as follows: 

&~U>° + £g • * : =-ISS™ ; ^=4{ dt ' e 4^'>' (29) 

The integrated form of the integrability condition, Eq. (|17p. then allows to express Xx> 
again through the averaged scalar curvature, Xp = Qkj) — Qz> — (R) v , and we obtain 
the sources corresponding to (|25[). however, with a curvature source that captures the 
deviations from a constant curvature model: 

*v - I \ - Qv - - gfcp/a|] ^ _ Q v [(H) v - 6k v /a 2 v ] 

In the following section we present new families of exact solutions of the effective 
Einstein equations. While models can be obtained on any chosen simply-connected 
spatial domain T>, most of the following solutions are global and the spatial domain 
T> is extended to the whole Riemannian manifold S, which we assume to be compact. 
(Note that regional solutions such as the example presented in Eq. (|27j) "appear" to 
only depend on the matter distribution inside T>, but this is not the case, since the 
initial data are to be constructed non-locally from the whole distribution in E; global 
solutions and their parameters therefore have a more robust status. 



3. Globally stationary effective universe models 

3.1. Globally static cosmos without a cosmological constant 

On the global scale we first require the effective scale-factor as to be constant on some 
time-interval, hence ds = as = and Eqs. (fT2*|) and (|13p may be written in the form: 

My 

Qs = 4ttG-4 - A ; (31) 
Vial 

(ft) s = 12ttG^ + 3A , (32) 

with the global kinematical backreaction Q s , the globally averaged 3-Ricci curvature 
(72) E , and the total restmass Ms contained in S. 

Let us now consider the case of a vanishing cosmological constant: A = 0. 
The averaged scalar curvature is, for a non-empty Universe, always positive, and the 
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balance condition (|3T|) replaces 0, while the condition replaces (jlj). Obviously, 
backreaction ()31|) and averaged scalar curvature (}3*2)) trivially satisfy the integrability 
condition (jlfijl . Thus, in view of (|25jl . the globally static inhomogeneous cosmos without 
a cosmological constant is characterized by the cosmic equation of state: 

= 3Q S = const. p e E ff = g^ s = . (33) 

3.2. Interlude: kinematical backreaction as a cosmological constant 

We first note that, apparently, in Eq. (JT2J) Qv plays the role of a positive cosmological 
constant; however, in Eq. (fT3*J) Q-p has opposite sign! That the kinematical backreaction 
term may take the role of the cosmological constant has been suggested in J7j and 
discussed in connection with the backreaction problem in j2Hj as well as in the recent 
discussion on Dark Energy and backreaction |H7| ITfij . The above remark shows that 
caution is in order with a direct identification. 

Of course, we may force the kinematical backreaction term to take exactly the role 
of the cosmological constant: neglecting A in the general equation (|T2j) . Q v may be 
regarded as a (in this equation possibly time-dependent) cosmological term. Then, 
from Eq. (fT3j) a constraint equation on Q v follows, if we identify all sources that 
imply a deviation from Friedmann's equation with the cosmological term; the correct 
requirement can be inferred from Eq. (fTTIJ) and reads: 

1 ^dt' Qv^-aUt') = Q v , (34) 



a 2 v J tl " " dt' 

which implies Qx> = Qv{U) — const, as the only possible solution. Physically, it appears 
contrived to freeze fluctuations to a constant value in an evolving model. However, it is 
noteworthy that Qx> indeed is required to be constant, if we force it to behave like A. 
Eq. (jl7)l then implies for the averaged curvature: 

= ~ 3Qv(U) • (35) 

The solution in the case kx> = has been noticed in |HZj too; however, freezing also the 
averaged curvature in an evolving model appears to be even more contrived. A different 
interpretation in the framework of an effective scalar field may imply a more meaningful 
interpretation of this solution (see Subsect. 13. 6|) . We finally write the cosmic equation 
of state for this special (regional) solution: 

P§f QviU) ~ ky/al 

gf s 8ttG {q) v (U)/al + Qv(U) - 3k v /a 2 v ' 1 ' 

which, for large \av\, approaches the equation of state p^ s = —gf s . 

3.3. Local instability versus global stability 

The problem of instability of the Einstein cosmos, as outlined in the introduction, 
must be thought of in two ways: i) a static homogeneous cosmos is unstable within 
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the class of Friedmann-Lemaitre cosmologies, since for a small change in density the 
balance between the density and the cosmological constant in Eq. Q is destroyed, 
leading to acceleration if AitGqh(U) > A, and to decceleration otherwise. Suppose we 
take a slightly smaller density, then Eq. (J3J) implies that H(t) > 0, and the model 
starts to expand; A (being constant in time) cannot dynamically compensate for this 
expansion, which itself decreases the density further, and so gives rise to the instability; 
ii) a homogeneous cosmos is unstable to inhomogeneous density perturbations. Such 
perturbations are amplified as a consequence of the attractive nature of the gravitational 
field tending to increase overdensities and to decrease underdensities. This is the 
content of local gravitational instability: inhomogeneities are amplified. The latter 
instability also applies to the other Friedmann-Lemaitre cosmologies. Both types of 
scalar instabilities (and also vectorial and tensorial perturbations) have been recently 
clarified and detailed by Barrow et al. [B] for the Einstein static universe model 
containing a perfect fluid, generalizing earlier work by Harrison [SU] and Gibbons |47j . 
[IE]. Further insight was added by Losic and Unruh [HI] who investigated the stability 
analysis to second order in scalar and metric fluctuations. A related analysis concerning 
the dynamical phase space as well as attractor or repellor properties of homogeneous 
solutions may be found in [33], [75] . 

Now, let us look at the same type of perturbations in the framework of the globally 
static, but inhomogeneous cosmos. Altering the density source would equally disturb 
the balance by virtue of Eq. (|3T|) . but it would not necessarily destroy it: the reason 
is that the kinematical backreaction term Qs indeed reacts back on this perturbation 
(justifying its name). For, unlike A, Q s can acquire a time-dependency and is then 
dynamically coupled to the averaged scalar curvature (TZ) E (in general through the 
integrability condition (|TB|l and in particular, for a stationary cosmos, through Eq. (|4"U|) 
below): as soon as the density is perturbed, taking for example a slightly lower average 
density as above, the (positive) curvature starts to evolve such that d t ( (7£) s V^ 3 ) < 0; 
in turn, d t ( QyYj) ) > 0, and, thus, fluctuations tend to decrease less rapid: as we shall 
see in an exact solution below, the coupling to the averaged scalar curvature can be 
strong so that decreases in proportion to the inverse volume similar to the density, 
not (as expected for fluctuations) in proportion to the square of the inverse volume. It 
does not matter if we consider homogeneous or inhomogeneous perturbations (as far as 
scalar perturbations are concerned), since the effective equations (fT2|) and (|T3|) govern 
both and are not narrowed to the class of homogeneous-isotropic solutions as in the 
case of the standard Einstein cosmos. In other words: a perturbed effective cosmology 
stays within the same class of cosmologies governed by the effective equations, since the 
latter are general. 

An interesting future task will be to analyze in detail - along the lines of a stability 
analysis of the homogeneous cosmos with symmetry [H], [M] - whether kinematical 
backreaction tends to stabilize perturbations of the dynamical balance between the 
globally averaged density source AnG {g) s and the global kinematical backreaction term 
<2s due to its coupling to the averaged 3-Ricci curvature (7^) s . In any case, we do not 
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expect that a generic cosmology would dynamically approach this state for any initial 
setting. In particular, a perturbed Friedmannian state may not approach a globally 
static state, rather the question is whether perturbations of an already established 
balance would destabilize this state. This calls for an investigation of perturbation 
theory of the global out-of-equilibrium state, rather than of a Friedmannian state. 



3.4- Globally stationary effective cosmologies 

Suppose that the Universe indeed is hovering around a non-accelerating state on the 
largest scales. Still, the effective static cosmos discussed above may appear as a quite 
rigid model. We may think that a more natural condition would be stationarity. Indeed, 
at first sight the balance condition, if attained, does not necessarily imply that the model 
is static. We may look for a wider class of models that balances the fluctuations and the 
averaged sources by introducing globally stationary effective cosmologies: the vanishing 
of the second time-derivative of the scale-factor would only imply ds = const. =: C, 
i.e., = as + C(t — ti), where the integration constant as is generically non-zero, 
e.g. the model may emerge [HHl, |10j from a globally static cosmos, as := 1, or from 
a 'Big-Bang', if as is set to zero. (In Friedmannian cosmology such an expansion law 
would correspond to a "curvature-dominated" model, since C 2 + k rs for negligible 
sources, hence resulting in a Hubble expansion that is determined by a constant negative 
curvature.) 

On the global scale and for a stationary cosmos, Eqs. (|12|l and (|13jl read: 



My 

Qs = 4ttG-4 - A ; (37) 
Vial 



(K) r = 12ttG^ + 3A - 6#| ; # s = — . (38) 



By inserting (J37)l into (|S5jl. we can evaluate the constant C by looking at the resulting 
equation at initial time; for the normalization as(ti) = 1 we get: 

6C 2 = 6A + 3Q E (*i) - (ft) E {U) . (39) 

We are now going to discuss the stationary cosmologies more explicitly by deriving an 
exact solution to the effective Einstein equations. (One can easily show from (J2J) that 
a Friedmannian cosmology does not allow for a stationary cosmos, since C 2 + k = Aa 2 
only allows for a = const.) 

First, note that the time-derivatives of the above equations deliver a dynamical 
coupling relation between and (TZ) s as a special case of the integrability condition 
(Hi): 

1 4C 3 

-d t Qy + -d t {n)^ = — . (40) 

In view of the conservation of the total restmass Ms, we infer directly from Eq. (}3*Tj) 
that Qs evolves as: 

& = _ a+ &M±a. (41) 
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For the same reason we infer from Eq. (J38|) that (7?.) E evolves as (inserting C 2 , Eq. (|39j0 : 
W ^3A + ( ^ (i -'- 3 2 fe(f - ) - 6A + 3fe(f - 3 )+3A . (42) 



The solution { (JSJ), } satisfies the integrability condition (fT7)J) : it provides the first 
example of a non-trivial solution to the effective Einstein equations, in which kinematical 
backreaction and averaged scalar curvature are dynamically coupled. (This solution, for 
A = has been first discussed in connection with the Dark Energy problem in [T5].) 

This solution stabilizes the (always positive-definite) combination 6A + 3Qs — 
(K) E , which evolves as (6A + 3Q E (^) - (H) s (ti))/°l = 6C 2 /4 = 6#|, so that 
for large \a^\ it approaches zero. This relation governs the "cross-talk" between 
kinematical backreaction, averaged scalar curvature and A. Although the stationary 
state approaches the condition needed for a static state, the time-dependency of the 
individual terms changes the global picture drastically. For, if we prepare an initial state 
in the vicinity of the globally static cosmos, the condition ()37|) is conserved in time, but 
the curvature evolves away from the condition (|38j) . which we shall explicitly discuss 
below for the subcase of a A— free stationary cosmos. 

To work with the stationary models has an advantage: we can employ the 
dimensionless characteristics (|2()jl . which are valid on the global scale, and so ease the 
discussion of observational results. Directly from the stationarity conditions ([37)1 and 
f!38|) we obtain for the dimensionless functionals: 

= -4fig + 2ttl and = 1 + 3fig - 3fi^ . (43) 

A straightforward calculation, employing the solution { (141}) . ()42j) }, provides their 
explicit solutions in terms of two parameters - chosen to be the mass density parameter 
at initial time and the parameter for the cosmological constant at initial time: 

ni. == n£(t.) = 2(gs ^ + A) ; ni := njfc) = A . (4 4) 

Inserting the constant C, Eq. (JBlJj) . we get: 

nl = ^- ng = _i^ + inX; ^ = i-^-|^4; «S = nX • (45) 

The cosmic equation of state for a globally stationary cosmos can be obtained by 
inserting (jSTJ) and (J3EJ) into (1231) : 



Pcff 1 

=: u'cff = -- 



1 - 3£l\ a| 



/ \ 1 



j 4A 



(46) 



3Q B (ti)-{K> s (ti) 

which is time-dependent for a non-vanishing cosmological constant and, although it 
approaches, for large |as|, the cosmic state p cS = —g e f[, the asymptotic model is in a 
stationary and not in a de Sitter phase, since the cosmological constant is assumed to 
share the global balance condition. Below we shall assign another role to A that will 
allow for a de Sitter phase too. For this purpose let us first drop A as a possible source 
sharing the stationarity condition. 
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3.5. Globally stationary cosmos without a cosmological constant 

Consider now again the subcase of a vanishing cosmological constant. The cosmic 
equation of state for a globally stationary cosmos with A = can be inferred from 
Eq. flUj): 

Pes = -3 Qes • (47) 

It is interesting to compare this condition with the investigation of backreaction in 
inhomogeneous cosmon fields by Christof Wetterich jHH] , in particular with the cosmon 
equation of state, which also points to a possibly interesting different interpretation of the 
backreaction term (cf. Subsect. I3.6|) . As already mentioned in the Introduction, in the 
case of a physical pressure source, the above equation of state allows for a static model 
without cosmological constant also in the homogeneous-isotropic case ([SI], p. 383). 
Recently, the cosmic equation of state in the present formalism has been calculated 
in the framework of the Tolman-Bondi solution [HH|, also confirming a strong coupling 
between averaged scalar curvature and kinematical backreaction. For related discussions 
see |2E|, US], EH- 

The solution for this A— free cosmos, as a subcase of the general solution { (JHJ), (JI2J) }, 

reads: 

y s = —3— ; Wy. = 72 + —3 — • ( 48 ) 

u s u s a s 

The total kinematical backreaction QsV-s, = 47rGMs is a conserved quantity in this case. 
This solution implies for the dimensionless functionals: 

nl = ^ ■ = ; fi| = i-^ ; n? = o, (49) 

as 4as 4as 

where this cosmos only depends on a single parameter, which we have chosen to be the 
restmass density parameter at the initial time £j, as in Eq. (|44|) . 

As already mentioned in the general case, the averaged curvature evolves strongly 
in this cosmology. In general, i.e. if the initial state is not the globally static state, it 
evolves away from the curvature of a static cosmos: 



2 



= 1 _ 6C 1 

12ttG(o) e 127rG(o) s a| 



1 - 



a s • (50) 



12ttG (q)x (U 

Stating this remark in terms of the solution (|4*2^) for A = 0, 

a S a S 

where the numerators in both terms are positive-definite, we find that the second term 
will dominate after some time, and the averaged curvature can change sign in the coarse 
of evolution. In view of (|4"Hj) . this will happen when fi^ = 1 — = 0, i.e. when the 
density parameter has dropped to the "critical" value = 4/3. 

A comparison of this cosmology with a Friedmannian cosmology in terms of the 
evolution of the single parameter Vt l m is performed in Appendix B. 
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3.6. Inflationary cosmogonies evolving into a globally static cosmos 

The inhomogeneous globally static cosmos and also stationary cosmologies were built 
on the assumption that an exact balance between averaged material mass and its 
kinematical fluctuations (not its density fluctuations) is established. The cosmo logical 
constant, if present, takes the role of a further player in an established globally balanced 
state: if kinematical fluctuations vanish, then the Einstein static model arises, if we have 
complete balance due to kinematical fluctuations, then A = 0. All intermediate states 
are conceivable. 

Now, instead of preparing such a non-accelerating state, we could further widen 
our model assumptions by allowing for a non-zero global acceleration. To illustrate 
this we shall choose the simple assumption that this acceleration is entirely due to the 
cosmological constant, and also that the corresponding expansion rate is also only due 
to A. In this case, the balance condition is established only among the averaged material 
mass and the kinematical backreaction as in the globally static model. We shall see that 
this "driven balance" (with the same cosmic equation of state as in the static model) 
is conserved and the fluctuating matter is subjected to a constant acceleration: the 
conditions (f3~Tj) and (|3^|) are assumed to hold (for A = 0), but the respective terms 
evolve in time according to the solution: 

; (e> ^<%M i (52) 

a s a E a s 

i.e., averaged scalar curvature and kinematical backreaction both obey conservation 
laws similar to the averaged density, and the "driven balance" is maintained. With our 
assumptions, a E (t) is given by the exponential solution of the flat de Sitter model. (Note 
that (72) E and Q E with (72) E = 3<2 E solve the integrability condition (fTB|) irrespective 
of the particular form of H s = as/a^.) 

The idea of inflation is related to matter creation. In this line a global stationarity 
condition singles out a state which assumes that, as soon as matter is created, also its 
kinematical fluctuations are large and are trying to establish the balance condition on 
the global scale. To address this question properly, the above model is too simple. A 
more ambitious model would first attempt to understand kinematical backreaction and 
averaged curvature by means of effective inhomogeneous scalar fields. We shall now 
discuss some aspects of such an attempt. 

In the absence of matter and extrinsic curvature fluctuations, Eqs. (JEJ), Cdj) and 
(fTSJl imply (we always consider the global scale here): 

3^= A ; 3(^V + <%M = A, (53) 
a E V a s/ 2a E 

leaving the exponentially expanding de Sitter cosmos as the only solution, when the 

initial scalar curvature is put to zero: H\ = a/A/3. As a next step, let us consider 

a cosmos without matter, but with intrinsic and extrinsic curvature fluctuations; the 
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effective equations (|T3|l and ()16)) then read: 

3^ = A+Q S ; 3f^V + ffk±gg = A . (54) 
as V a s/ 2 

As in ()26|) the above equations can be recast into standard zero-curvature Friedmann 

equations with new effective sources, now implying an interpretation in terms of an 

effective scalar field: 

11 11 

Q$ ■= Qcs = --^^Qs - rr— ^ (ft) s ; P<s> ■= Pes = --^^Qs + (ft) E • (55) 

lb7rG 16/tG 1d7tG 487tG 

In this matter-free case the above equations suggest an interpretation of the geometrical 
degrees of freedom in the extrinsic curvature fluctuations as the kinetic energy part of 
a scalar field: Qs obeys a stiff equation of state; also, depending on the sign of Qs, we 
may have "phantom energy" (see also e.g. [72]). An important remark here: the 
simplified case of vanishing averaged scalar curvature implies through the integrability 
condition (jTB} that <2s = Qz(ti)a^. e , and the effective sources decay in proportion 
to the square of the inverse volume. As the previous considerations about globally 
stationary cosmologies have shown, the average curvature plays a significant role for 
the maintainance of a large Qs and therefore, the presence of curvature is crucial and 
should not be neglected in this picture. 

If we now take matter sources into account, the cosmic equation of state changes 
and it is here, where detailed models of matter creation would have to be analyzed to 
understand whether and if, physically, such a creation process would entail, e.g. strong 
global expansion fluctuations, favoring an inhomogeneous state at the time when A no 
longer rules the expansion. If Qs is interpreted as a kinetic energy term of a scalar 
field (we may call it a morphon), then its very presence would also shape the extrinsic 
curvature of created matter fluctuations; its role in the equations above shows that a 
possible balance condition in the process of conversion of kinetic energy into matter 
would be <2x = 47rG (f?) s , as required for a globally static or stationary cosmos. 

3. 7. Global far- from- equilibrium states 

Expressing a fluctuation-dominated cosmos in a thermodynamic language, the globally 
stationary state is in an out-of-equilibrium state compared with the Friedmannian 
"equilibrium state" (in the sense defined below). In [53J an entropy measure has 
been advanced that we can employ to characterize both states: looking at the non- 
commutativity relation (jHJ) for \I/ = g on the global scale, Hosoya et al. |HS1 found 
that the source of non-commutativity is given by the production of relative information 
entropy, defined as to measure the deviations from the average mass density due to the 
development of inhomogeneities: 

(d t g)x-d t (g)x = - ^ , (56) 

with the, for positive-definite density, positive-definite Kullback-Leibler functional 

S{q\\(q)x}--= I Jd 3 X glnj^- . (57) 
Jt, \Q)y, 
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This measure vanishes for Friedmannian cosmologies ("zero structure") and so defines 
the notion "equilibrium state" introduced above. It attains some positive time- 
dependent value otherwise. The source in (JHH)) shows that relative entropy production 
and volume evolution are competing: commutation can be reached, if the volume 
expansion is faster than the production of information contained within the same 
volume. 

In j^S] the following conjecture was advanced: 

The relative information entropy of a dust matter model S{g\\ (q)^} is, for 
sufficiently large times, globally (i.e. averaged over the whole compact manifold Sj 
an increasing function of time. 

This conjecture can be proven for linearized scalar perturbations on a Friedmannian 
background (the growing-mode solution of the linear theory of gravitational instability 
implies d t S > and S is, in general, time-convex, i.e. d 2 /d 2 S > 0). However, in 
a generic non-perturbative situation it may not hold. Below we shall approach this 
question for the case of a globally stationary cosmos. 

Let us first consider the general situation. We calculate the second time-derivative 
of (|57jl and first obtain 

'8q5B) + (gh(A6) 2 , (5* 



\ We 

where, for any scalar field 5^ := \I/ — denotes the deviation from the 

global average value, and A\l/ := a/ (<5\I/) 2 ) s the global amplitude of its fluctuations. 
Raychaudhuri's equation specifies the deviations 58 for the time-evolution of the 
expansion rate, 

S6 = -AnG5g - -5{6 2 ) - 25{a 2 ) , 5A = , (59) 
3 

which, together with the commutation rule yields: 

§- = AnG(Ag) 2 + (^) s (A6) 2 + J (5g56 2 )^ + 2 (5g5a 2 )^ . (60) 

For our purpose the above equation may be recast by explicitly writing out the 
last second terms, and expressing the expansion fluctuation amplitude in terms of 
Qs = f (A6>) 2 - 2 (a 2 ) E (cf. Eq. JUJ)). We obtain the general equation: 



4" = B i + (Q)j 



1 



(61) 



3 

with B 1 := 4ttG(Aq) 2 + | (gO 2 )^ + 2 (ga 2 }^ > 0. 

From this equation we infer that, for a globally static universe model with = 
47rG {g) s and (0) s = 0, the Kullback-Leibler functional is always time-convex. 

For the case of a stationary universe model (with or without A) we can obtain 
sufficient conditions for time-convexity as follows. Directly from (|61|) we conclude that 
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Qs — | = Qs — 3 iff, should then be positive and, employing the dimensionless 
characteristics (J2Uj) on the global scale, that the condition 

^2 < ~\ (62) 

is sufficient. Inserting the stationarity condition = AttG (f?) s — A into (pTTj) we instead 
need the condition AnG (g)^ — A — | (9)^ > which, expressed through the average 
characteristics (}2T)j) on the global scale, reads: 

inS-nJ > i . (63) 

For a A— free stationary cosmos we conclude that a sufficient condition for time- 
convexity of the Kullback-Leibler functional is ift > 2, which implies a sufficiently 
positive global average curvature through Eq. (f4*3~)): 

n| = i - lift =► n| < -1 . (64) 

Given the above results, we have three possibilities for the evolution of the 
entropy functional (|H7j) in a globally static cosmology (and also in a globally stationary 
cosmology, if the conditions above are met), depending on the sign of its first time- 
derivative, the relative entropy production rate (Eq. note that for the special case 
of a globally static cosmos df (g) E — 0): 

firstly, if S is positive at initial time ti, then the functional ()57|) is always growing; 
secondly, S becomes positive at some time t > ti, in which case the functional (|57j) 
is always growing thereafter; thirdly, S never becomes positive, i.e. the functional 
(JH7J) approaches some constant Sq > from above. This latter case would violate the 
Entropy Conjecture advanced in |HS1- In the first two cases, looking back to Eq. (fHHjl. the 
static model leaves the volume unchanged, while information is created. This enhances 
the source of non-commutativity, which may be interpreted as a signature of a non- 
Friedmannian state. 

The globally stationary cosmos does not necessarily imply a time-convex evolution 
of the relative information entropy. In this context it is interesting to think of a 
floating equilibrium realized in open biological systems (since gravity is long-ranged, 
a gravitational system is also not isolated). In this analogy a stationary far-from- 
equilibrium may be characterized, according to Prigogine [73 , by a minimum of the 
entropy production rate, i.e. for our inhomogeneity measure (157)1 we would require the 
condition: S = 0. Whether such a saturation is possible within a dust cosmology is 
doubtful in view of a lack of physics like pressure gradients. However, the above analogy 
prevails. 

3.8. Discussion: the global picture 

An investigation of models in the framework of the effective equations while still 
including a cosmological constant shows that the inflationary picture of cosmogonical 
scenarios in the Early Universe has not to be abandonned. The cosmology can still 
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undergo several phase transitions and inflationary phases. Like in the early days after 
Einstein's suggestion of the static cosmos, it is a subject of controversy related to the 
initial conditions of the universe model: either the singular 'Big Bang' initial state, or 
the static initial state eventually evolving into an expanding universe model is preferred 
(see the discussion in [8J. In the context of inhomogeneous inflationary scenarios, 
we are entitled to single out an averaged model that is globally stationary, i.e. we 
assume that inflation may end and, during a phase of damped oscillations, exit into a 
globally stationary phase, instead of a homogeneous Friedmannian phase. If such an 
inhomogeneous "initial condition" at the end of inflation is prepared, then e.g. a globally 
static inhomogeneous cosmos certainly is in a stronger position than his homogeneous 
predecessor with regard to its more general properties, especially if its stability could 
be proven. Similar ideas have actually been already discussed within chaotic inflation, 
e.g., (S3 ESI, being mirrored here within a classical cosmos. 

Justifying such an inhomogeneous initial state physically is another issue. If we 
follow Einstein's thoughts leading to the static cosmos, the relevant arguments were 
not empirical, but rather philosophical. These thoughts are far from conceiving a 
'fitting model' to observational facts. The cosmological principle that elevates the 
Friedmannian models to a matter of principles is not always the underlying reason 
to employ homogeneous-isotropic cosmologies. As the stationary but inhomogeneous 
cosmology shows, there is no local isotropy and, hence, the cosmological principle must 
be replaced by another principle of "global prejudice" (in the constructive way of the 
word) like for example: "the Universe as a whole cannot create momentum out of 
itself" (a "non-Munchhausen" principle + ). The formulation of a sound principle is 
itself a considerable task; the tentative balance condition that was helpful to construct 
exact examples, should be replaced by a cosmic virial theorem averaged on the global 
scale, which in turn calls for a generalization of the dust matter model. Also, global 
topological constraints should be important. The way how to deal with a cosmological 
constant, e.g. the question whether we should include it in a virialized state [HZj, and the 
more fundamental approach of an effective scalar field dynamics in the Early Universe is 
another issue of concern in the formulation of a global principle. For example, a global 
net acceleration can be a result of initial conditions, and in fact this is the case in generic 
models of Friedmannian cosmology. 

While the globally static cosmos features properties similar to his homogeneous 
predecessor, e.g. the global density and the global curvature are constant in time, 
the globally stationary cosmos evolves very differently. We have shown that, although 
the relation among kinematical backreaction and averaged curvature tends to the 
corresponding relation in the globally static cosmos, the time-dependency of the 
individual terms may affect the global properties strongly. The stationary state tends 
to the static state only in the sense that, e.g. in the case of an expanding cosmos, the 
rate of expansion slows down, but the steady increase of the scale factor allows for a 

+ Baron v. Munchhausen was trying to pull himself out of the pond without external help. 
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global change of the sign of curvature. As Eq. (|51|) shows, an initially positive averaged 
scalar curvature would decrease, and eventually would become negative. This may not 
necessarily be regarded as a signature of a global topology change, as a corresponding 
sign change in a Friedmannian model would suggest; the averaged scalar curvature 
is only a weak descriptor for the topology in the general case, and information on 
the sectional curvatures is required for definite conclusions on this issue (see jl] Ch.9; 
0). However, if a global change of topology arises, we could imagine that the cosmos 
started with a finite volume of a closed spaceform, then the averaged curvature of 
the inhomogeneous space may become negative and a globally hyperbolic curvature 
arises that would suggest a finite-volume compact spaceform with (now no longer 
simply-connected) sections of negative sectional curvatures and, generically, a "horned 
topology" (for topology-related issues see [EH, jHHl, [H] 

3. 9. Discussion: the regional picture 

The implications that are furnished by the effective static model featuring regional 
expansion and contraction are qualitatively very different from a cosmogony with global 
expansion: firstly, the former is more "violent" than the latter, since fluctuations grow 
exponentially rather than in a moderate power-law pace; secondly, an expanding region 
must be counterbalanced by a contracting region. The second property implies that, 
depending on estimates of the total size of the cosmos, we should be close to seeing 
a large "blue excess" in galaxy number counts, e.g. at the borders of the currently 
drawn Sloan survey. The microwave sky would be redshifted in expanding domains 
alternating with a blue-shifted sky in contracting domains, while we would have to 
claim an underdensity in our regional Hubble volume [HU EH] , |EH] , counterbalanced 
by overdense regions beyond the horizon; we would be surrounded (in space and time 
direction) by high-density walls. This spatially and temporally alternating scenario 
implies a large paradigmatic change compared with the standard model of cosmology. 

Exponential fluctuations imply the likely situation that, e.g. the global curvature 
cannot be seen on the scale of our Hubble volume. Not unlikely, since we are seeing 
a large Hubble expansion, the regional curvature is (relatively) more negative, i.e. in 
a special case seemingly zero or negative. As already noted, a globally static cosmos 
would necessarily call for a replacement and, regionally for an obvious refinement of the 
cosmological principle including the question, whether and how close our observers have 
to be at the center of such a regional "Hubble bubble" jHS], j^H], 0. The scale of this 
"reduced curvature region" likely exceeds scales that have been discussed in connection 
with peculiar-velocity catalogues that are statistically affected by boundary conditions 
PT] , A strong constraint is the isotropic microwave sky, a question related to 

assigning the whole dipol to our proper motion [2H] • 

The globally stationary cosmology evolves very differently and, hence, also its 
regional properties would be different. Since already the global curvature can approach 
0, the expanding stationary cosmos may show similar regional fluctuations compared 
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with a Friedmannian model, however, in a much narrower range of times. The differences 
to a Friedmannian model can be clearly seen in the detailed comparison given in 
Appendix B. 

We shall not embark here further into speculations about implications of a globally 
static or stationary cosmos. It is clear that such a drastic change of views stimulates 
thoughts, but detailed calculations needed for proper statements do not exist yet. 
Rather, we shall put a currently held discussion about the Dark Energy problem into 
perspective, which turns out to lie at the heart of the balance condition required for 
globally stationary cosmologies. 

4. The Dark Energy problem 

Let us first recall the contemporary view on the global universe model. The standard 
model of cosmology idealizes the matter distribution in the Universe to be homogeneous 
and isotropic, neglecting structure. Observations point to an averaged matter content 
of at most 30 percent (including dark matter) and almost flat space sections jSOj- The 
missing gap to fill the global cosmic triangle (jA.2|) is modeled by a (spatially constant) 
cosmological term in the simplest case [H] (the sum of the parameters has to be equal 
to 1). This is the concordance model (roughly 0.3 + 0.0 + 0.7 = 1), which can be 
fitted to a large set of orthogonal observational data. (There are, however, other voices 
[T3*] . [T2*].) Modeling this dark energy gap with a cosmological constant results in an 
accelerated phase*. The 'coincidence' that dark energy starts to dominate exactly when 
also structure enters the non-linear regime suggests that there could be a physical 
relation between the effect of structure on the average expansion (known as backreaction 
effect) and the dark energy gap found in the standard model, thus providing a natural 
solution to this coincidence problem. (It is clear here that sub-horizon fluctuations are 
made responsible.) 

From Eq. (|12jl the condition for an accelerating patch T> of the Universe directly 
follows: 

Q v > AnG{g) v . (65) 

This regional condition is weaker than the requirement of global acceleration, since 
it accounts for the regional nature of our observations. It has been recently claimed 
[H7j that the condition (jHSJ) could be satisfied within our regional Hubble volume, 
hence providing a smart explanation of the Dark Energy problem [SB], In JH] 

the following argument has been advanced. The considerations in this paper and the 
"classical" explanation of the Dark Energy problem in terms of backreaction effects are 
intimately related through the condition (|65p. If this condition would hold, then also the 
cosmologies presented in this work would attain the status of physically viable models 

* The only direct support for acceleration comes from high-redshift supernovae observations |78j : 
however, the 'fitting models' for the interpolation between low-and high-redshift observations and the 
measured distances do not take inhomogeneities into account, see |38| for a discussion and references. 
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on the global scale, since the physical basis of these models is exactly the possibility that 
kinematical expansion fluctuations and averaged sources are of the same order. A crucial 
support for the condition to survive in time could also be given in ^B] employing 
a particular example of the averaged Einstein equations, listed as globally stationary 
cosmos without A in this paper. This example of an exact solution indeed points to a 
strong coupling between kinematical backreaction and averaged scalar curvature: the 
rate of decay of Qs is in proportion to (f?) E and can therefore be of the same order as 
4tcG (g)j, today. 

Caveats 

There is, however, a large body of opponents including myself who do not think that the 
condition (|65J1 can be attained within the standard model of cosmology. However, this 
opposition rests on the picture that fluctuations grow through gravitational instability 
of a homogeneous universe model. If we allow for the very different picture of 
inhomogeneous initial conditions, then the above outlined explanation can be a natural 
consequence. In the standard picture of cosmological structure formation from Cold 
Dark Matter initial conditions there are essentially three caveats that have to be 
overcome, which we are going to put into perspective now in some more detail. 

Caveat 1 : explicit calculations of kinematical backreaction 

The above summarized suggestion of solving the Dark Energy problem has a long 
history; there have been many attempts to calculate backreaction effects following the 
advent of the averaging problem initiated by George Ellis |35] (see [H], [131 EH HDl EE] , 
[77] , [TH E5J EH] , |BB]; [II], [SJJ and refs. therein). The new input into the discussion 
concerns the quantitative importance of backreaction summarized by the "bold claim" 
(JHSJ)- Earlier efforts to support this claim point to a negative answer. Let us therefore 
look into the details of explicit calculations of the effect. 

Cosmologists are mostly thinking in Newtonian terms when modeling structures 
(e.g. by N-body simulations which are constructed such that their spatial average 
evolves as in the standard model). Computation of backreaction in the Newtonian 
framework is possible in great detail [23]; but there is a drawback (often giving rise 
to misunderstanding in the literature): the relevance of kinematical backreaction for 
the global evolution cannot be estimated within Newtonian cosmology, simply because 
one can in general prove that an averaged Newtonian cosmology has zero global 
kinematical backreaction [22] Q.v%- Nevertheless, as non-perturbative calculations in 

jj The reason is that the kinematical backreaction term Qx> can be written as a total divergence on 
Euclidean space sections, and a Newtonian cosmology, in order to be uniquely defined, has to impose 
periodic boundary conditions on the inhomogeneities [22122; hence, by Gauss' theorem, backreaction 
terms are just boundary terms and they vanish for an empty boundary! Putting this fact into a 
constructive perspective, this non-trivial property of Qx> proves that the architecture of current N- 
body simulations is physically correct. 
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[2U] show, backreaction has indirect impact on the evolution of the standard cosmological 
parameters on (relatively large) regional scales. However, the dimensionless backreaction 
parameter (|27Ijl f2g := — on large expanding regions is quantitatively negligible 
compared with the effective density parameter := 3 y^ G a ¥x'i , i-e. the condition 
rewritten in terms of dimensionless characteristics, 

-n%>^, (66) 

(implying for the effective deceleration parameter (|22|) < 0) is not fulfilled on large 
scales in these calculations. In fact, VLq is not only small, but positive, i.e., with regard to 
([14)1. shear fluctuations dominate over expansion fluctuations. In [23] it was achieved to 
fulfill this condition on a patch of 100 Mpc (for h = 0.5) by saying that, on these smaller 
scales, we may have different expectation values for velocity fluctuations (the kinematical 
components entering Q v in Newtonian theory) and for the density fluctuations (they 
are theoretically independent and, on 100 Mpc, it may well happen that the expectation 
value for Qz> is 1-cr off the expectation value for the density fluctuations, i.e. we 
could prepare an 'untypical' sample). Although this situation is likely to happen (the 
probability distribution of Q-o is non-Gaussian), we cannot apply this argument to very 
large samples. 

To give another example (out of many that are not quoted here), the answer 
attempted in the work by Russ et al. [ZZj was based on a second-order perturbative 
calculation in General Relativity. However, initial and boundary conditions were 
assumed such that, already from their assumptions, the global backreaction term, 
according to our definition, has to vanish; see the discussion in JH]). Although it is 
obvious from the above remarks about Newtonian calculations that the global value 
can only be determined in the framework of General Relativity, the quoted relativistic 
calculation was "designed" after a Newtonian model. This also reveals inherent 
difficulties for treating the problem with sufficient generality, if standard assumptions 
like periodic boundary conditions are adopted. 

Given these remarks we could say that Caveat 1 is due to Newtonian or quasi- 
Newtonian work, or (since calculations have been mostly done perturbatively), is 
generally due to the fact that backreaction effects need a non-perturbative and 
background-free calculation (see jHH] for suggestions of non-perturbative variables in 
General Relativity). Until more general calculations have been done, we may say that 
the above quoted results cannot be extrapolated. However, basically a non-perturbative 
calculation in General Relativity has already been done: we have pointed out that 
Eq. ()19|) for the scale-factor is identical to the corresponding Newtonian equation. The 
non-perturbative Newtonian model investigated in j2Sj was based on integrating this 
exact equation for a specific fluctuation source term, the (perturbative) 'Zel'dovich 
approximation'. Employing instead the relativistic version of this approximation as 
given by Kasai (SEj does not lead to differences concerning the results either. Therefore, 
such a relativistic calculation for a© yields the same result in terms of Q-p. The difference 
is that the latter is linked to the scalar curvature (and this may give rise to different 
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interpretations as will be examined in a forthcoming work [24 ). But, the quantitative 
conclusions given in j2H| can to some extent be extrapolated. This remark also points to 
the necessity of considering initial conditions other than those assumed in [2Bj in order 
to meet the condition (JH3J). 

Caveat 2: observational constraints 

The effective cosmological parameters defined in (|2T)J) can be considered to provide a 
fair representation of the values which also an observer would measure in a sufficiently 
shallow survey region T> (the light-cone effect is not taken into account). We may 
therefore discuss estimates of those parameters in comparison with observed values. 
Note, however, that the interpretation of observations is mostly done by employing 
a standard Friedmannian cosmology as a 'fitting model' and therefore, geometrical 
inhomogeneities (that are hidden in the definition of the spatial averages in the 
Riemannian volume element) are ignored; below we shall argue why we should not 
ignore them. 

On the grounds that the parameters ()2())1 can be "observed", we examine the 
condition which itself would "only" require a contribution of |f2g| > 0,075 for 

the value = 0, 3 to the cosmic quartet. Assuming that this condition holds on some 
large domain V, which we may take to be as large as our observable Universe, then 
Hamilton's constraint in the form (}2T|) also implies: 

nj + = i - n£-ng > 1 - , (67) 

showing that we need a substantial negative curvature (positive fl^) on the domain 
T>, if we put the cosmological constant to zero. (To reconcile this condition with 
a non-vanishing cosmological constant would need an even larger value of fl^ than 
that suggested by the concordance model.) The fact that a large value of kinematical 
backreaction entails a substantial average Ricci curvature has also been stressed by 
Rasanen [7B] . 

The condition (|57jl seems to contradict the widely agreed expectation that the 
curvature should be very small, the concordance model assumes an exactly zero scalar 
curvature. However, the Friedmannian curvature parameter is only directly related 
to the averaged scalar curvature in the homogeneous case. If we would - mistakenly - 
estimate the averaged scalar curvature through = —k v /{a^H'^) ) then we would have 
to estimate the (regional) time-integral of Qv in the relation fl^ + Qf = 1 — VL^ 1 — VLq N 
(c/. Eq. ([23])), and thus our estimate would depend on the dynamical evolution of the 
kinematical backreaction parameter and not just on its value. (In the (global) example 
presented in Appendix B we would have Q% = 1 — Q l m = const., while the true curvature 
functional is strongly time-dependent!) 

If we ask whether the kinematical backreaction Qx> is observable, the answer within 
our setup above is yes: on the observable domain D, Qd is built from invariants 
of the peculiar-velocity gradient in a Newtonian model. Ignoring again geometrical 
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fluctuations (and with them the fact that in a relativistic setting Qx> cannot be 
represented through invariants of a gradient, which is derived from a vector field) good 
high-resolution maps of peculiar-velocities could in principle determine the value of 
backreaction. Existing catalogues are, however, too small and they would only return 
the cosmic variance around the assumed Friedmannian background in a likely untypical 
patch of the Universe. 

On the other hand we know several observational facts that could place constraints 
on the value of kinematical variables jSEj- "Global" bounds on Qx>, where T> is of 
the order of the CMB (Cosmic Microwave Background) scale, can be inferred from 
work of Maartens et al. [65] . [66] : observational upper bounds were given in terms 
of covariant and gauge invariant quantities, e.g. for an upper bound on the shear 
fluctuations (o- 2 ) v we use the limit cr 2 /9 2 < 16a 2 10~ 10 , where both a* and 9* have 
to be interpreted as averaged values in our setting, i.e. cr* := {o~) v and := 3Ht>; the 
parameter a is determined by observations and depends on scale; for large-scale CMB 
observations we set a ~ 1. In the case when shear fluctuations would dominate, we may 
assume {e) 2 v = {6 2 ) v . In that case, Qx> reduces to —1(o 2 ) v (is therefore negative); the 
corresponding dimensionless parameter obeys the bound 



< g^ = «c < 2,4 a lO" 9 for (9) 2 V = (0 2 ) v . (68) 

This bound is comfortably satisfied in the calculations quoted above, but it is not 
what we need for explaining Dark Energy: the opposite situation must be examined, 
namely, we assume shear fluctuations to be negligible and search bounds for a positive 

Qv^U(e 2 ) v -(e} 2 v ). 

On the Hubble length L H = 2c/ H « QGpc/h Maartens and his collegues [HE] 
obtain an upper bound on the density fluctuations that we may convert into expansion 
fluctuations using, for simplicity, the proportionality relation stemming from the 
standard linear theory of gravitational instability on an Einstein-de Sitter reference 
background with density qh and Hubble expansion parameter H (suppressing the 
decaying mode): 



(Q-{Q)v)l . 



) 2 
vlv 



Qh K 2 H 2 ' (69) 

where the initial volume Vi = 4/37r(L#/2) 3 accounts for the dimension and the scale. 
With these assumptions and the upper bound on the density fluctuations quoted in [66J, 
the backreaction parameter is bound from below: 

°>i§ = ^>-(^ + s) 2 * 2ur8 for {a% ^' (70) 

i.e., for VL^ = 0, 3 on the Hubble-scale we would have VLq > —8, 9 a 2 1CT 7 . 

However, a precise relation of the bounds on gauge-invariant variables given by 
Maartens & collegues with the relevant variables discussed here requires more efforts and 
is the subject of a forthcoming work. Notwithstanding, the above order-of-magnitude 
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estimates demonstrate that observational constraints will also provide obstacles to 
advocate the condition 

The loophole here: the strong coupling of density fluctuations to expansion 
fluctuations, as given by a direct proportionality for small perturbations of a 
Friedmannian model, Eq. (|69j). was strongly used to place constraints on expansion 
fluctuations. In an out-of-equilibrium state these two types of fluctuations (being 
independent variables) are decoupled and an observational constraint must be based 
directly on an upper bound on expansion fluctuations. 

Caveat 3: geometrical fluctuations, curvature backreaction and the choice of foliation 

It is a reasonable question to ask whether a cosmological model is accelerating or not 
given the framework of the kinematically averaged Einstein equations, provided we 
consider the chosen foliation of spacetime appropriate. Caveat 3 consists of the fact 
that those averages are performed on an inhomogeneous hypersurface. This is all right 
as long as model and observations are to be compared, since observations themselves are 
carried out in an inhomogeneous geometry. The problem we speak about here comes 
into the fore by 'fitting' a Friedmannian cosmology to observational data, and this is 
the standard procedure. The question whether a cosmological constant is needed to 
explain the observations is due to this fitting process. Consequently, there is a need to 
compare averages of the dynamical model to averages on a homogeneous geometry, which 
is a difficult task because the averaging of tensors is not straightforward jHSl EH EE] 
There is obviously an interpretation problem in the standard model (having Euclidean 
geometry in the case of the concordance model): averages on inhomogeneous spaces 
have not the same values as averages on homogeneous spaces. Riemannian volumes 
may substantially differ from Euclidean volumes and are not small perturbations per se 
(a Riemannian ball has 7r 2 /6 less volume than an Euclidean ball with the same geodesic 
radius [21] )• It is possible to define a "geometrical renormalization" flowing averages 
on inhomogeneous spaces into Friedmannian averages ^Hj. The consequences for the 
regional cosmological parameters are summarized in [2U] implying further "backreaction 
effects" such as Ricci curvature backreaction that lead to a "dressing" of cosmological 
parameters by smoothed-out geometrical inhomogeneities. These effects may provide a 
further key for explaining Dark Matter and Dark Energy. 

The answer to the question whether these effects are significant needs detailed 
realistic models for "volume roughening" of spatial slices in the Universe, and robust 
estimators for intrinsic curvature fluctuations. Since this problem has been largely 
overlooked when considering averages in cosmology, there are not much results. A naive 
swiss-cheese model by just glueing Riemannian balls in place of Euclidean balls in the 
slice yields an effect of 67% reducing the necessary dark energy roughly from 70% to 
50% [21]. A mismatch of similar magnitude has been reported by Hellaby [22] using 
volume-maching as suggested by Ellis & Stoeger [JT] (compare [HI])- However, the 
quantitative significance of these effects on cosmological scales have yet to be explored. 
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Related to the averaging of both extrinsic and intrinsic inhomogeneities 
Zalaletdinov [HD] proposed a macroscopic description of gravitation based on a covariant 
spacetime averaging procedure. The geometry of the macroscopic spacetime follows 
from averaging Cartan's structure equations, leading to a definition of correlation 
tensors. Macroscopic field equations (averaged Einstein equations) can be derived in 
this framework. Within this approach Coley et al. |27j have recently shown that for 
a spatially homogeneous and isotropic macroscopic spacetime, the correlation tensor is 
of the form of a spatial curvature term. In this context it was also speculated that the 
extra correlation terms might help to stabilize a globally static universe model. 

As for any averaging procedure the choice of spatial slices on which one considers 
average quantities is crucial. In ^H] this problem has been identified and, in the 
framework of smoothing the geometry with the help of the Ricci flow, one has strong 
means to approach the slicing problem, e.g. by aiming at "minimizing" artificial gauge 
effects [1^]. Recently, Ishibashi and Wald [55] gave particular examples that show why 
caution is in order concerning a straightforward observational interpretation of averaged 
kinematical quantities on given spatial slices. Here, a clear advantage of volume- 
averaging scalar quantities is that the resulting averages of 4-covariant variables and 
their governing equations can be transformed to a different choice of slicing, which is a 
feature of the covariant fluid gauge [37] underlying the averaged equations in the present 
paper; examples that illustrate this have been given in 

5. Concluding Remarks 

In this work we have contrasted Friedmannian cosmology to particular choices of 
inhomogeneous cosmologies. Two very different points of view concerning the 
comparison of an idealized Friedmannian cosmos and an averaged inhomogeneous 
cosmos are conceivable: 

Firstly, an averaged inhomogeneous model could be a small perturbation of a 
Friedmannian model in the sense that the kinematical backreaction on the global scale, 
Qs, is small and, hence, the kinematics of both models measured by the evolution of 
the effective scale factor a-s are comparable. Still, as shown in detail within Newtonian 
cosmology in [2H], this must not imply that the parameters of a Friedmannian model 
evolve as in the inhomogeneous model; on the contrary, even for negligible backreaction, 
there are significant differences. We may place such a universe model into a near- 
equilibrium state as measured, e.g. by the information-theoretical measure (j5T|) that 
vanishes for a Friedmannian cosmology. 

Secondly, and this point of view was exploited in the present work, we may conceive 
a universe model in a far- from- equilibrium state, characterized by strong averaged 
expansion fluctuations on the global scale. Such a model is no longer a perturbation of 
a Friedmannian model. We have shown that these highly-inhomogeneous models may 
then evolve in the vicinity of the balance condition <2 S = A-nG (g)^ (disregarding the 
cosmological constant). 
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Which of these two points of view is realized depends on the choice of the initial 
state of the Universe at the exit epoch of an eventual inflationary phase in the Early 
Universe. We have argued why a possible solution of the Dark Energy problem through 
kinematical backreaction effects is possible from the second point of view, but unlikely 
from the first point of view. Given such a global state we have shown that strong 
kinematical fluctuations are conserved, and the state is maintained by a strong coupling 
of backreaction to the averaged Ricci curvature. This has been exemplified by an exact 
solution of the averaged Einstein equations for a globally stationary cosmos. We may 
say that large kinematical fluctuations are maintained on the cost of averaged scalar 
curvature in this solution. 

We have considered three families of cosmologies that belong to the out-of- 
equilibrium state of vanishing global acceleration: the first family covers globally 
stationary cosmologies without a cosmological constant; as a subcase, a globally 
static cosmos is possible that revives the original Einstein cosmos concerning its 
global properties, however, having very different properties on regional scales. The 
second family covers models that obey the stationarity condition of vanishing global 
acceleration, but include a cosmological constant sharing the balance condition. In 
this case we have also given the exact solution. Finally, the third family rests on the 
stationarity conditions of the first family, but the global acceleration is non-zero due 
to the cosmological constant alone. Fluctuations and average characteristics in such 
A— driven cosmologies have very similar properties to the previous cosmologies, however, 
with the presence of a global exponential expansion. 

Although, formally, these new families of cosmologies obey simple equations and 
are in this respect very close to Friedmannian cosmologies, their dynamical properties 
must be regarded as significantly different. Investigating details of the new models 
must be considered as an endeavour of high magnitude, especially if one is interested 
in the regional properties. One important reason for this difficulty lies in the scale- 
dependence of the average characteristics. For example, the Einstein static model has 
globally and locally a constant curvature and a constant density, while the globally 
static, but inhomogeneous cosmos features strong matter and curvature fluctuations on 
regional scales. The latter model enjoys the freedom that the averages can be scale- 
dependent, e.g., the Einstein static model has (like all Friedmannian cosmologies) a 
scale-independent density, while the globally static cosmos allows for a decay of the 
average density by going to larger scales. It is possible that fluctuations, encoded in 
Qz>, on a regional domain V are large relative to the averaged density 4nG (q) v on that 
domain, but Qy, could be very small on the global scale |E| in absolute terms, and still 
of the order of AttG (f?) E . 

These new models need detailed investigation before any conclusive result 
relevant to observations can be obtained. In light of the present investigation 
three lines of research appear to be fruitful strategies: firstly, investigating global 
stability properties of the presented cosmologies, i.e. studying perturbations on non- 
Friedmannian "backgrounds" including studies of gravitational instability on regional 
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domains; secondly, investigating the relativistic generalization of [22! and comparing 
standard with globally inhomogeneous initial conditions [24J and, thirdly, investigating 
inhomogeneous inflationary scenarios and the properties of fluctuations at the exit 
epoch, i.e. studying the role of kinematical backreaction in the Early Universe. 
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Appendix A: the size of a globally static cosmos 

In order to get an idea about typical numbers characterizing a globally static cosmos, 
we are going to discuss some numerical estimates. We already imply that there exist 
regional fluctuations of this global state; we take the classical Einstein cosmos as the 
global model and infer our estimates from regionally observed properties. 

Dividing Friedmann's differential equation (0) by H 2 (here, the static case is 
excluded), we obtain Hamilton's constraint for the homogeneous-isotropic case, written 
in the iconized form of a cosmic triangle [7j: 

Vt m + + Vt k = I ; Vt m := ^}p H ' ^ A := 3^2 ' ^ k := a 2}j2 ' (A-2) 

Confining ourselves now to our regional Hubble volume, we write f2 m ( Hubble— scale) — : Ct m 
and express H through h := H/lOOKm/sMpc. From (jSJ) with k = +1 and (JUJ) we first 
obtain: 

a E = i. ^ '- = —ZGMe = — CISchwarzschild • (A-3) 
^/AiiGqe TT TT 

Extrapolating the value of the matter density g H = ^^a m H 2 = qe, we find for the 
Einstein radius in space units: 

r E = ca E = -^=rVV3 1(T 8 ^ , (A.4) 

which, for a m ~ 0, 3 and h ~ 0, 6 yields r E ~ 7, 5 Gpc. For comparison: if a m = 1 and 
/i ^ 0, 46 [12] we would get r E ~ 5, 3 Gpc, for a m w 0,3 and h w 0,46, r E ~ 9,7 Gpc. 
Defining a Hubble volume roughly by Vh = ^-(c/H) 3 (where we took the horizon 
radius of a Friedmannian model, r# := c/H w 3Gpc/h, as radius of an Euclidean 
sphere randomly placed within the Einstein cosmos), we get for the first set of values 
above Vh ~ 174, 5 Gpc 3 , and with Ve = 2ir 2 r E ~ 8319 Gpc 3 , we conclude that a 
fluctuating static cosmos with the radius of the homogeneous Einstein cosmos would 
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roughly contain 50 Hubble volumes. Thus, the total volume V# is comfortably large to 
allow for significant regional fluctuations on the Hubble scale. 

The regional curvature that we would observe on the Hubble scale will not be the 
global curvature of the Einstein cosmos. However, if we assume for simplicity that the 
curvature is non-fluctuating, but that we have a fluctuating expansion being globally 
zero and regionally given by h - as will be possible in the general model discussed in 
the text), we could estimate the global Friedmannian curvature parameter from regional 
observations: setting flk( Hubble— scale) — '. Olfo — —k/af^H 2 , we obtain with the first set of 
values above a k = 10~ 16 s 2 Gpc 2 jm 2 « —0,44. 

Thus, if we would 'fit' a Friedmannian model within the regional Hubble volume, 
we would obtain from (jA.2|) (0, 3; —0, 44; 1, 14) for values of the regional cosmic triangle 
(a m ; ak, «a), i-e. straight application of this (quickly cooked) model would be in trouble 
compared with the values of the so-called concordance model as a 'best-fit' Friedmannian 
cosmology in comparison with observations, (0, 3; ±0; 0, 7). 

However, these estimates are naive in a variety of ways: firstly, the averaged 
curvature is not given by the Friedmannian curvature parameter, if kinematical 
fluctuations are present; secondly, a fluctuating cosmos implies fluctuations in scalar 
curvature too, i.e. it is easily conceivable that the magnitude of the regional curvature 
could be smaller than the above estimate. Moreover, as we have learned, if expansion 
fluctuations are assumed large, then we can as well dismiss the cosmological constant 
altogether; thirdly, a further feature of globally static inhomogeneous cosmological 
models is that they would allow for a significantly larger size compared to the values 
quoted above, since the average density may likely be scale-dependent and may decrease 
towards the global scale. 



Appendix B: evolution of ft— parameters in the globally stationary cosmos 
compared with their evolution in Friedmannian cosmologies for A= 

At first glance, the reader may think that a globally stationary cosmos, introduced 
in Sect. 3, looks very similar to a Friedmannian cosmos, since the scale factor has a 
powerlaw form, as oc t, and the global expansion H^, oc t~ l , but in fact both cosmologies 
are drastically different in nature. This can be nicely illustrated with the evolution of 
the density parameter, which is the only free parameter in both cosmologies, if A = 0. 

For this comparison we employ Eq. ([19)1 on the global scale, and evaluate this 
equation with the solution of the stationary cosmology Qs = Qn{U)a^?: 

H 2 (l - fl£)a| + A; s = | £ dt' Q s a' E a s = ~Qx(U) (J- - lj , (B.2) 

with fcs = (Qs(ti) + (7£) s (U))/6 = H 2 {VL l m — 1), and the density parameter fl^ in the 
general model (|2T)J) . From the definition of the latter we have a second equation: 



2 



H 

fi m = ^mjrnj ■ ( B - 3 ) 
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Inserting (jB.3|) into (|B.2|) . a simple calculation provides: 

nE _ l^n /r> A\ 

m ^(l-K) + [l-^(l-«)]a s ' 1 • ] 

with k = 1 for the stationary inhomogeneous model with non-zero kinematical 
backreaction, Qs(ij) = ^Kn-^-h an d ^ = for the Friedmannian models. This minor 
formal difference amounts to a substantially different evolution; the reader may insert 
numbers to compare the following solutions: 

fl m (Friedmann) = — , ~ tt^ r J fl (stationary) = — — . (B-5) 

fi m + (1 - WJax " a s 

Introducing the redshift z in the Friedmannian model, as = and accordingly a 

(formal) effective redshift z efr in the inhomogeneous model, we can write the above 
solutions in terms of the density parameter today, Cl l m := fl^(z = 0), as follows: 

Qi (I _i_ z ) 

^(Friedmann) = ^ -. ' J ^(stationary) = fi' (1 + 2 cS ) . (B.6) 

1 + ft^z 

The scalar curvature in a Friedmannian cosmology, fl*. = 1 — Q m cannot change sign, 
but the averaged scalar curvature in the globally stationary cosmos, 
(see Eq. can. (Note, that the "Newtonian" curvature functional, introduced in 

Eq. (jUJ), stays constant in the stationary solution: = 1 — f2JJ) 

While the density parameter in a Friedmannian model stays smaller (larger) than 
1, if it is initially smaller (larger) than 1, and while the curvature parameter stays 
smaller or larger than determined by initial data, the corresponding parameters (and 
the additional backreaction parameter) in the stationary model can "communicate with 
each other" and can "freely operate" (of course, subjected to the Hamiltonian constraint, 
+ &n + = 1) an d the stationarity conditions) in the range of values: 

4 2 1 

+ oo > tt^ > +- ; -oo < < +- ; -oo < fig < -- . (B.7) 

ODD 

(These conditions can be easily derived from the fact that 3Qs + (^)s = > 0.) 

The infinities are reached at the "Big Crunch" of a contracting stationary model with 
C < 0, in which case as — > in a finite time. 

It should be emphasized that these are the global values of the cosmological 
parameters; on regional scales very different values may be measured, whereas in the 
standard model of cosmology global and regional values have to be equal. 



